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I The classical Bondi-Hoyle formula gives the accretion rate onto a point par- 

ticle of a gas with a uniform density and velocity. However, the Bondi-Hoyle 
problem considers only gas with no net vorticity, while in a real astrophysical 
^ I situation accreting gas invariably has at least a small amount of vorticity. We 

I therefore consider the related case of accretion of gas with constant vorticity, for 

the cases of both small and large vorticity. We confirm the findings of earlier two 
dimensional simulations that even a small amount of vorticity can substantially 
change both the accretion rate and the morphology of the gas flow lines. We 
show that in three dimensions the resulting flow field is non-axisymmetric and 
time dependent. The reduction in accretion rate is due to an accumulation of 
circulation near the accreting particle. Using a combination of simulations and 
analytic treatment, we provide an approximate formula for the accretion rate of 
gas onto a point particle as a function of the vorticity of the surrounding gas. 
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1. Introduction 

Accretion of a background gas by a point-like particle is a ubiquitous phenomenon in 
astrophysics, occurring on size scales ranging from black holes in galactic nuclei to individual 
stars or compact objects accreting the winds of their companions. Bondi, Hoyle, and Lyttle- 
ton (Hoyle and Lyttleton 1939, 1940a,b,c; Bondi 1952) treated this situation by considering 
the ideal problem of a massive point particle within a uniform gas moving at a constant 
speed with respect to the particle. For a polytropic gas, the problem is largely characterized 
by a single dimensionless parameter, the Mach number Ai, which is defined as the ratio of 
the gas-particle velocity to the sound speed of the uniform gas. (The ratio of specific heats 
7 also affects the accretion, but its effects are only order unity.) The Bondi-Hoyle formula 
gives the approximate accretion rate as a function of A4, which, in a version updated by 
Ruffert (1994) and Ruffert and Arnett (1994), is 

Mbh = 47rpooG'2M2c-3 

Here, A is a number of order unity that depends on M. and on the equation of state. For 
M. — (Bondi accretion) in an isothermal medium, A — exp(1.5)/4 « 1.1. 

The Bondi-Hoyle formula finds wide application in astrophysics. Our first motivation 
for exploring this topic is star formation. The competitive accretion picture of star forma- 
tion posits that protostars form with initially low masses, typically ^ O.IM©, and grow in 
mass by accreting unbound gas from the molecular clump in which they form (Bonnell et 
al. 1997; Bonnell, Bate, and Zinnecker 1998; Bonnell et al. 2001a,b; Klessen and Burkert 
2000, 2001). The competition for gas is invoked to explain the initial mass function (IMF). 
This is a process of Bondi-Hoyle accretion, albeit in a turbulent medium, and one ought 
to be able to use a Bondi-Hoyle-like formula to estimate the rate at which the proposed 
"seed" protostars would accrete. A second significant area of application is in black holes, 
particularly the supermassive black holes (SMBHs) at the centers of galaxies. Numerous 
authors have used the Bondi formula to estimate the rate at which SMBHs accrete, and 
found that it substantially overestimates the accretion rate (Di Matteo, CariUi, and Fabian 
2001). Proga and Begelman (2003) have argued that rotation may explain the difference 
between accretion rates estimated from observations and those predicted by the Bondi rate. 
A third area of application is in accretion rates onto compact objects. Perna et al. (2003) 
suggest that a combination of rotation and magnetic fields leads to lower accretion rates onto 
isolated neutron stars than one would estimate using the Bondi formula. They argue that 
this lowered accretion rate accounts for the unexpectedly small number of isolated neutron 
stars that have been detected. 
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Previous authors have considered the problem of hydrodynamic accretion with non-zero 
vorticity, both analytically (Sparke and Shu 1980; Sparke 1982; Abramowicz and Zurek 
1981) and numerically (Fryxell and Taam 1988; Ruffert 1997, 1999; Igumenshchev and 
Abramowicz 2000; Igumenshchev, Abramowicz, and Narayan 2000; Proga and Begelman 
2003). Abramowicz and Zurek (1981) analytically treated accretion of gas with non-zero 
angular momentum onto a black hole. They found that for very low angular momentum the 
flow pattern is similar to Bondi accretion; the flow is quasi-spherical, it becomes sonic at 
a radius much larger than the Schwarzschild radius, and the accretion rate is comparable 
to the Bondi rate. For high angular momentum, the gas forms a disk and the sonic radius 
in the equatorial plane becomes only a factor of a few times the Schwarzschild radius. The 
accretion rate is substantially below the Bondi rate. Proga and Begelman (2003) performed 
two dimensional axisymmetric simulations of the accretion of slowly rotating gas onto a black 
hole. They find three regimes of accretion. For sufficiently low total angular momentum, the 
fiow is Bondi-like, similar to the pattern Abramowicz and Zurek (1981) predicted. For inter- 
mediate total angular momentum, the highest angular momentum gas that cannot accrete 
forms a dense torus about the black hole. The torus blocks off accretion except through a 
thin funnel along the poles, reducing the accretion rate to a value roughly independent of the 
total angular momentum. Finally, for even higher total angular momentum, a centrifugal 
barrier, coupled with the dense torus, further reduces the accretion rate. 

Our goal in this paper is to unify this work by describing the accretion process using a 
one-parameter family similar to the Bondi-Hoyle-Lyttleton treatment, and to give an anal- 
ogous formula for the accretion rate as a function of this parameter. We do this using a 
combination of three-dimensional simulations and analytic treatment. In these two impor- 
tant respects, our work is more complete than previous work. In §2 we give an analytic 
analysis of the problem, and present our approximate formula. In §3 we describe the com- 
putational methodology we use in our simulations, and in §4 wc compare the results of our 
simulations to the analytic theory we put forth. We summarize our conclusions in §5. 



2. Analytic Treatment 

We consider a point particle accreting in a medium with non-zero vorticity. The simplest 
case of this, which we treat here, is a medium with constant density and vorticity far from the 
particle, and where the gas has no net velocity of relative to the particle. Thus, we consider 
a particle at the origin of our coordinate system, surrounded by a gas whose velocity far 
from the particle is 

y 

Voo = Ci;*Coo — X, (2) 
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where tb = GM/ is the standard Bondi radius, M is the mass of the particle, Coo is the 
sound speed of the gas at infinity, and y is the y coordinate. The factor uj^ is a dimensionless 
number which we call the vorticity parameter, since the vorticity of the velocity field is 
—uj^Coo/'rBZi independent of position. 

The significance of a;* is clear: — 1 corresponds to the case where gas with an impact 
parameter y = t-q is traveling at the sound speed at infinity. The Keplerian velocity at a 
distance of tb from the central object is 

/ X fGM 

vkKTb) = \ = Coo. (3) 

V 

Thus, for cu^ — 1, gas whose impact parameter is tb is traveling at the Keplerian velocity 
for that impact parameter. Since vk oc r~^/^, while for our initial velocity field v^o oc y, gas 
with smaller impact parameters is sub-Keplerian, while gas with larger impact parameters 
is super-Keplerian. For a;* < 1, the transition from sub- to super-Keplerian occurs at 
impact parameters larger than r^; for a;* > 1, the transition occurs at impact parameters 
smaller than tb- Because the Bondi radius is the distance at which gas comes within the 
gravitational sphere of infiuence of the central particle, provided its kinetic energy at infinity 
is small, we should expect a transition in behavior at a;* f« 1. 

We note that previous authors (e.g. Abramowicz and Zurek (1981); Proga and Begelman 
(2003)) have usually parameterized their flows in terms of speciflc angular momentum rather 
than vorticity, as we do here. The choice of parameterization is somewhat a matter of taste. 
Specific angular momentum has the advantage that it yields a circularization radius for the 
fiow that depends only on the specific angular momentum, simplifying the analysis. On the 
other hand, the choice of vorticity enables us to invoke the Kelvin Circulation Theorem in our 
analysis, which we find helpful. In addition, since specific angular momentum is a function 
of accretor position, and gas farther than tb from an accretor is unaware of its presence, 
arranging a flow with constant speciflc angular momentum requires an unlikely coincidence. 
It is more natural to characterize a flow by its velocity gradient, hence its vorticity, which 
is independent of accretor position. Regardless of the choice of parameter, however, it is 
possible to make an approximate translation of the formulae we derive based on a;* to ones 
in terms of the specific angular momentum. For the initial conditions we are considering, 
the specific angular momentum of a gas parcel with impact parameter y is 

1/2 

loo = yVoo = i^*Cs ■ (4) 

Tb 

In the artificial case of a fiow with constant specific angular momentum rather than constant 
vorticity, the equivalent value of a;* should be determined roughly by the value of the vorticity 
at the Bondi radius. Thus, as an approximate conversion we suggest a;* = Ioo/{csT^b)- 
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Before moving on to analyze the behavior of accretion as a function of a;*, we consider 
briefly the question of the gas equation of state. We choose to restrict ourselves to consid- 
ering an isothermal gas, and we denote the constant sound speed as c^. This is physically 
appropriate for the star formation context: radiative energy transport in the dense molecu- 
lar gas from which stars form is capable of eliminating any temperature variations on times 
small compared to the mechanical time scales in the gas (Vazquez-Semadeni et al. 2000). 
In the context of SMBHs and compact objects, the thermodynamics is more complicated, 
as the gas can range from isothermal to adiabatic depending on the efficiency of radiative 
energy transport. However, based on prior simulation work with Bondi accretion (Ruffert 
1994; Ruffert and Arnett 1994), and accretion with shear (Ruffert 1997), we do not beheve 
that changes in the equation of state will change either the flow pattern or the accretion rate 
by more than a factor of order unity. 



Abramowicz and Zurek (1981) performed the flrst detailed analysis of the cases of small 
and very small a;* (which we define below), that we shall review and extend. We can 
distinguish two regimes that occur for c^;* <^ 1, depending on the physical dimensions of the 
accreting object. In the absence of viscosity or other mechanisms of angular momentum 
transport, a parcel of gas has constant specific angular momentum l^o as it falls towards the 
accretor. The orbit of the gas parcel changes from predominantly radial to predominantly 
tangential at a characteristic circularization radius, where the parcel's velocity is roughly the 
Kepler velocity. The specific angular momentum of a parcel of gas in Keplerian rotation at 
a distance r from the central object is 



and the circularization radius is therefore determined by the condition l^^, — iKit^drc)- Equat- 
ing (4) and (5), a gas parcel at initial impact parameter y has a circularization radius 



If the circularization radius for a given parcel of gas is small compared to the radius of the 
accretor, r^, then that gas will reach the surface of the accretor traveling in a predominantly 
radial orbit. Its angular momentum, and its vorticity, will therefore be deposited in the 
central object and will not otherwise affect the flow pattern, leading to an accretion rate 
approximately equal to the Bondi rate. 



2.1. Very Small a;* 



lK{r) — rvxir) = vGMr, 



(5) 




(6) 
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In the case of small vorticity, where the velocity of the gas at infinity is small compared 
to the Keplerian velocity, gas begins to feel the influence of the accretor at a characteristic 
radius of r-Q. We can therefore distinguish our two regimes based on the relative sizes of 
^^€(^"6) and ta- If ?"circ(^B) < ''"Ai then we are in the case analogous to that of a small 
circularization radius in a flow with constant specific angular momentum: the gas that 
reaches the accretor will be on primarily radial orbits, and the accretion rate should be the 
Bondi rate. The condition for this to occur is 



We term the case tu* < c<;crit "very small a;^,." To get a sense of values of cjcrit in typical 
cases, consider some of the applications we discussed in the introduction. In a star forming 
region, a typical accretor mass is M ~ Mq and a typical sound speed is ~ 0.2 km s^^, 
giving a Bondi radius of re ~ 0.1 pc. A typical protostellar radius is ta ~ 10^^ cm, so 
Wcrit ~ 10"^ For the SMBH at the center of the Milky Way, M a; 3.7 x lO^M© for the 
standard distance estimate (Ghez 2004). The temperature is not constant and depends on 
the details of radiative processes, but a reasonable estimate for gas far from the black hole 
is ~ 10^ K (Meha 1994), giving Coo ~ few x 10^ cm s~^ and tb ~ 0.1 pc. The accretor radius 
is the Schwarzschild radius ta = 2GM/c^ ~ 10^^ cm, so again a;crit ~ 10~^. 

Note that this condition uj^ ^ a;crit for Bondi-like accretion is roughly equivalent to the 
conditions of Abramowicz and Zurek (1981) and Proga and Begelman (2003) on the specific 
angular momentum when the accreting object is a non-rotating black hole. In that case, 
the accretor radius is the Schwarzschild radius, ta = 2GM/c^, so a;crit = •\/2cs/c, and the 
specific angular momentum of gas with impact parameter y = tb is 



In their analytical treatment, Abramowicz and Zurck (1981) estimate that the critical specific 
angular momentum for Bondi-like accretion l^o = 2rAc; however, the exact critical value will 
depend on the details of how the accretion fiow joins onto the central object. 



If the vorticity is larger than ci^crit, but still small compared to unity, then gas circularizes 
before it reaches the accretor. The accretor is effectively a point particle, and cannot absorb 
angular momentum or vorticity from the flow. In this case, Proga and Begelman (2003) have 
found in their simulations that gas with too much angular momentum to accrete builds up 
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^oo(^b) = t^critC^rB = V^rAC. 



(8) 



2.2. Small 
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into a torus around the accretor. This thick torus inhibits accretion by blocking off a large 
part of the solid angle through which streamlines could otherwise reach the accretor. 

We can explain the growth of the torus by making a useful analogy to magnetohydro- 
dynamics. The flow at inflnity contains a constant vorticity. In regions of the flow that are 
non- viscous, the Kelvin Circulation Theorem requires that the vorticity contained in an area 
element of the fluid is constant. This is analogous to flux-freezing in a magnetic medium. 
The accrctor's gravity drags lines of vorticity in toward the accreting particle. Near the 
particle, where viscosity may become significant (for example due to magneto-rotational in- 
stability in the accretion disk), gas can slip through the vortex lines and reach the accreting 
object. However, as long as the region of significant viscosity is confined near the accretor, 
the vortex lines cannot escape. As accretion drags gas inward, vortex lines build up near the 
accretor, just as magnetic flux would build up in a magnetic flow. 

The effect of this vorticity build-up is to increase the rotational velocity in the gas. 
Consider a circular curve C of radius rc in the xy plane, centered at the origin. The 
circulation around curve C is the line integral of the velocity. 



where oj is the vector vorticity at any point in the flow and A is the area inside curve C. As 
a reference state, consider a velocity field with v = Vqo everywhere. The accretion process 
increases vorticity within the curve relative to this reference state, so it must also increase 
r and hence the mean rotational velocity around C. Because the vorticity parameter 
uj^ < 1, we know that initially < vk as long as rc < t-q. However, this analysis shows 
that increases with time, so eventually the rotational velocity will become comparable to 
the Keplerian velocity and a rotationally supported disk must form. The Keplerian velocity 
is supersonic for all radii smaller than tb, and the infiow onto the disk is not purely radial, 
so the disk will likely develop internal shocks and a time-dependent structure. 

The process of vorticity increase near the accretor cannot continue indefinitely. Eventu- 
ally, inward transport of vorticity by accreting matter must be balanced by an outward flow 
of vorticity in gas leaving the vicinity of the accretor and escaping to inflnity. The radius at 
which gas becomes unbound from the accretor and can start escaping is tb, and we therefore 
expect that the disk will extend out to ~ tb in radius. The scale height of a thin Keplerian 
disk is 




(9) 



and by Stokes' theorem 




(10) 
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where Q is the angular velocity. Since the disk radius is ~ tb, we expect if ~ tb- Thus, the 
scale height is comparable to the radius and we have, rather than a disk, a thick torus. The 
torus may be somewhat thicker for an equation of state with ratio of specific heats 7 greater 
than unity, because in that case the sound speed will increase inside the Bondi radius due to 
compressional heating. However, the compression should be is only order unity at r ~ Tb, 
so the effect of compression on the torus scale height should also be order unity. 

The circulation trapped within the torus should be 

r ~ 2'KrBVK{rB) = ^ttCsTb- (12) 

In our reference state with v = Voo everywhere, the circulation is 

To = 7rr|a; = irCgrBCO*, (13) 

which is clearly much smaller than the equilibrium value of F for c<j* ^ 1. The time required 
for the circulation to build up to F is determined by the rate at which mass falling onto 
the accreting object carries vorticity inwards. The gas within the Bondi radius reaches the 
central object, or the growing disk, on the Bondi time scale, ^b = ^b/cs. Thus, the flow 
replaces the gas within the Bondi radius on this time scale. Each set of "replacement" gas 
carries the same vorticity as the gas it is replacing, and therefore increases the circulation 
by Fq. As a result, the time required for the accumulating vorticity to reach its equilibrium 
value is approximately 

F 2 

^eq ~ TT^B = — ^B- (14) 
i ^* 

This estimate only applies for a;* ^ 1, when accumulation of circulation is the key factor in 
setting the final accretion rate. 

We can roughly estimate the amount by which the accretion rate will be reduced by 
the presence of the torus once the flow reaches equilibrium. The torus reduces accretion 
by blocking off solid angle through which material could otherwise accrete. If the torus 
extends out to tb, and its height is roughly tb as well, it blocks all angles in the range 
n / A < 9 < Sn / A. If the reduction in accretion rate is proportional to the fraction of solid 
angle blocked, then the accretion rate for this case should be reduced to ~ 30% of the Bondi 
rate, Mb, given by equation (1) for = 0. 

It is important to note that the equilibrium state we describe here does not depend on 
the value of a;*, as long as o^crit < tu* <C 1. The equilibrium radius of the torus depends 
only on the maximum vorticity that can accumulate near the accretor before the velocity 
becomes comparable to the Kepler velocity and halts infall in the disk. Thus, the equilibrium 
radius of the torus is ~ tb independent of cu^, and the scale height and accretion rate are 
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also if ~ tb and M ~ 0.3Mb, respectively. Different values of a;* only mean that vorticity 
builds up more or less slowly, and thus that there is a longer or shorter initial transient. 



2.3. Large a;* 

The final case to consider is a;* ^ 1. In this regime, the velocity of the gas at the Bondi 
radius is comparable to or larger than the Keplerian velocity. To derive an approximate 
accretion rate, we adopt an ansatz from Bondi-Hoyle accretion. 

Regardless of its velocity, gas at infinity initially has a positive total energy. Its potential 
energy is zero and it has non-zero internal energy. In order to become bound and accrete 
onto the point mass, the gas must lose its energy by shocking. If gas does shock, it is likely 
to accrete. Since the flow at infinity is laminar, shocks will only occur for those streamlines 
that are bent significantly by the gravity of the accreting particle. To determine if a straight 
streamline is susceptible to bending, we can compare the velocity of gas on that streamline 
(adding both bulk and thermal velocity) to the maximum value of the escape velocity along 
that streamline. If the peak escape velocity is larger, then the streamline will bend and is 
likely to terminate in a shock. Suppose the gas at infinity is traveling in the x direction, and 
consider a streamline that begins at xyz coordinates {±oo, y,z). The streamline will bend 
significantly and pass through a shock if 

2GM . .9 9 

>v^{y,zf + cl. (15) 



In the Bondi-Hoyle case, Voo — M.Cs independent of y and z. Thus, (15) reduces to 

2GM 



< ,2(1 + ^2) = 2rBH, (16) 

where we have used the standard definition for the Bondi-Hoyle radius. 

If we make the assumption that all parcels of gas traveling on streamlines that bend 
significantly go through shocks, and all gas that shocks loses enough energy to become bound 
and eventually reach the accretor, then the rate at which mass accretes onto the particle 
is equal to the rate at which gas from infinity begins traveling on streamlines that meet 
condition (15). We can therefore define an area A in the plane at infinity that satisfies (15) 
and approximate the accretion rate by 



M 



I P oo^/ V oo{y, + cl dydz, (17) 

J A 
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where poo is the density of the gas at infinity and we have again added the bulk and thermal 
velocities of the gas in quadrature. For the Bondi-Hoyle case, A is simply a circle of radius 
of radius 2rBH and v^^, is a constant, so the integral is trivial. Evaluating it gives 



Mbh ~ 47rr|HPooC.Vl + ^ 4^^^, 



which is just the Bondi-Hoyle formula (equation 1) up to factors of order unity. 

Because this ansatz allows us to reproduce the Bondi-Hoyle formula, we will apply it to 
predict the case of accretion with vorticity. In this case, — oj^CsU/r-Q, so (15) becomes 

2GM f CsV 2 2 

>[u:,-] y^ + cl (19) 



With the convenient change of variable y' = y/r-Q and z' = z/vb, this becomes 

2>{y" + z"y/\l+u;ly"). (20) 
We can solve (20) for z' to find the function Z'{y',u;*) that defines the boundary: 

We refer to region of the yz plane bounded by ±Z'{y', a;*) as A{uj^). The maximum value of 
y' on the boundary of A{u!^) occurs for z' — 0, and is 

2/3 



y'm^M = '—Tfs • (22) 



3V3 (qu;, + 73 + 81^) - 32/3 

Figure 1 shows the boundary of ^4(0;*) as a function of a;*. 
The approximate accretion rate is 

M{cu,) ^ rlpooCs [ ^/l + Z^dy'dz' (23) 

= AnrlpooCsfM, (24) 
where we have defined the numerical factor 

/(a;,) ^ ^[ ^l+^dy'dz' (25) 

VT+My^ / dz' dy' (26) 



- / \ U ^ (27) 
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The integral is straightforward to evaluate numerically, and we plot the result in Figure 2. 
We can also obtain an approximation for /(a;=„) in the limit a;* — > oo. We show in Appendix 
A that in the limit a;* ^ 1, we can approximate /(a;*) by 

/(^*) ~ ^ In (16^,) . (28) 
otto;* 

We also show this approximation in Figure 2. 

We must make one further modification to our estimated accretion rate. In the limit 
^* 0, f{w*) — > 1 and (24) reduces to approximately Mb. Thus, this approximation 
smoothly interpolates between the cases a;* -C 1 and a;* ^ 1. However, we have already seen 
that even in the case a;* <^ 1 the accretion rate can be substantially below the Bondi rate, 
since the accumulation of circulation leads to the formation of a thick torus that blocks 
streamlines from reaching the accretor. The same phenomenon should happen with a;* ^ 1. 
Circulation will still accumulate near the accretor. The shape of the torus should be the 
same as in the a;* -C 1 case, because that is set just by the physics of disks. We therefore 
reduce our estimated accretion rate (24) by a constant factor so that it matches our estimate 
of 0.3Mb for the case cUcrit < cu* <^ 1. The required factor is 0.34; it is not exactly 0.3 because 
(24) for a;* ^ 1 differs from the Bondi rate by a small factor. 

With this modification in place, (24) should give a good approximation of the accretion 
rate for all a;* > a;crit- Our final estimate for the accretion rate by a flow with vorticity a;* 
at infinity is 

M(.,).4.p„(^.|-i;(l^f/^ ■ (29) 

[ 0.34 /(a;*) : a;* > a;crit 

In §5 we give a more accurate approximation that is calibrated by our simulations. 



3. Computational Methodology 

To test the theory presented in the previous section, we ran a series of simulations. 
In this section, wc describe the simulation methodology, and in the next we compare the 
simulation results to our theory. 
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3.1. Code 

The calculations in this paper use our three dimensional adaptive mesh refinement 
(AMR) code to solve the Euler equations of compressible gas dynamics 

|^ + V-(pv) = (30) 

^(pv) + V-(pvv) = -VP-pV(t> (31) 

^(pe) + V-[(pe + P)v] = pvV0, (32) 

where p is the density, v is the vector velocity, P is the thermal pressure (equal to pc^ 
since we adopt an isothermal equation of state), e is the total non-gravitational energy per 
unit mass, and (j) is the gravitational potential. The code solves these equations using a 
conservative high-order Godunov scheme with an optimized approximate Riemann solver 
(Toro 1997). The algorithm is second-order accurate in both space and time for smooth 
flows, and it provides robust treatment of shocks and discontinuities. Although the code 
is capable of solving the Poisson equation for the gravitational field (j) of the gas based on 
the density distribution, in this work we neglect the gas self-gravity and include only the 
gravitational force of the central object. The potential is therefore 

, GM , , 

= , 33 

r 

where M is the mass of the central object and r is the distance from it. We do not adopt 
the Paczynski and Wiita (1980) potential frequently used for simulations involving black 
holes because we are do not wish to hmit ourselves to the case of black holes, and because 
even in the black hole case we are interested in the regime where the Schwarzschild radius is 
extremely small on scales of the grid, so that general relativistic effects are not important. 

Our code operates within the AMR framework (Berger and Ohger 1984; Berger and 
CoUela 1989; Bell et al. 1994), and is described in detail in Truelove et al. (1998) and Klein 
(1999). We discretize the problem domain onto a base, coarse level, denoted level 0. We 
dynamically create finer levels, numbered 1, 2, ... n, recursively nested within one another. 
To take a time step, one advances level through a single time step Ato, and then advances 
each subsequent level for the same amount of time. Each level has its own time step, and in 
general At;_|_i < At;, so after advancing level we must advance level 1 through several steps 
of size Ail, until it has advanced a total time Ato as well. In all the simulations we present 
in this paper, we chose cell spacings such that At; = 2Ai;+i, and thus we take two time steps 
on level 1 for each time step on level 0. (We find that refining by factors of two gives better 
accuracy in the solution than using a larger refinement factor.) After each level time step. 
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we apply a synchronization procedure to guarantee conservation of mass, momentum, and 
energy across the boundary between levels and 1. However, each time wc advance level 1 
through time Ati, we must advance level 2 through two steps of size At2, and so forth to 
the finest level present. 

On the finest level of refinement, we represent the central, accreting object with a sink 
particle (Krumholz, McKee, and Klein 2004). One noteworthy feature of our sink particle, 
which we shall discuss further in §3.3, is that the sink particle does not accrete any angular 
momentum from the gas in the computational grid. This is in contrast the sink methodologies 
used by previous authors (e.g. RufTcrt (1997); Proga and Bcgclman (2003)), where the sink 
could accrete angular momentum as well as mass. Our approach is appropriate for cu* ^ cUcrit, 
which is the case on which we focus. 



3.2. Initial and Boundary Conditions 

For each run, we place a sink particle at the origin of an initially uniform gas. The gas 
has an initial velocity field 

Vo = LO^Cs—X, (34) 

where we use different values of cj* in different runs. Table 1 summarizes the values of UJ^, 
that we simulated. Our computational domain extends from — lOOrs to lOOre in the x and z 
directions. In the y direction, we also used a range of — lOOre to lOOre for smaller values of 
a;*. For larger values of a;*, though, the large velocities that occur at large values of y produce 
very small time steps that make the computation prohibitively expensive. We therefore use 
a smaller domain in the y direction, as we indicate in Table 1; in every case, however, we 
chose our domain to extend to values of y such that i'o(z/max) ^ i'escape(z/max)- 

We used infiow/outfiow boundary conditions in the x direction and symmetry in the y 
and z directions. However, for all our runs the boundary is sufficiently far away from the 
central object that, within the duration of the run, no sound waves can propagate from the 
central object to the boundary and back. Our lowest resolution runs have Axmin < 
which we use for smaller values of a;*. For larger values of a;*, we use Axmin < rB/160. 
We set up our adaptive grids such that, at any radius r, the local grid resolution is Ax < 
max(r/20, Axmin). 
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3.3. The Sink Particle Method and Very Small Vorticity 

As Table 1 indicates, we vary uo^ from 10^^ to 10^'^, thereby thoroughly exploring the 
small uo^ and large uo^ regimes. However, none of our work explores the very small u;* 
regime. This is for two reasons, one physical and one technical. In order to explore the very 
small uj^ case, one must use either an unphysically small value of uj^ or an unphysically large 
accretor radius, leading to an unrealistically large a;crit- As we have shown in §2.1, for typical 
astrophysical situations in which accretion with vorticity is important, a;crit ~ 10~^. This is 
a truly tiny amount of vorticity, corresponding to a shear of one thousandth of the sound 
speed at one Bondi radius. In the case of a protostar in a molecular clump, this corresponds 
to a velocity gradient of no more than ~ 10 cm s~^ over a distance of ~ 0.1 pc; for the 
galactic center, it corresponds to a gradient of ~ 10'^ cm s^^ over a distance of ~ 0.1 pc. It 
is difficult to see how to produce such an irrotational flow field, and thus the very small a;* 
case is not relevant for the applications with which we are concerned. 

The technical reason we do not treat the very small cj* case is that our sink particle 
method is constructed to exactly conserve total angular momentum during the accretion 
process. (For details, see Krumholz, McKee, and Klein (2004).) Transfer of vorticity from 
the flow to the accretor is the distinguishing characteristic of flow with uj* < cUcrit- However, 
since our sink particle does not change the angular momentum of the flow fleld, it actually 
increases vorticity by removing mass while leaving angular momentum. Our code does 
dissipate vorticity via numerical viscosity, as we discuss in more detail in §4.4 on convergence 
testing. However, this effect is small except in the inner few zones. On balance, we consider 
this approach more realistic than the standard technique of allowing any angular momentum 
that enters a chosen accretion region to accrete. In reality, a small accretor should absorb 
negligible angular momentum from the flow on scales much larger than the accretor radius. In 
real accretion disks, viscosity acts to transfer mass inwards and angular momentum outwards. 
This is exactly the approximation we adopt in our sink particle method: mass travels inwards, 
angular momentum does not. While ideally one would follow the flow down to the true 
physical surface of the accretor, this is computationally infeasible in more than one dimension 
even with AMR. Our sink particle method provides a more realistic approximation of the 
behavior of an accretion disk than would allowing angular momentum to accrete. 



- 15 - 



4. Simulation Results 

4.1. Density and Velocity Fields 

In each of the simulations, the flow went through a transient and then settled down into 
a quasi-equilibrium state. The time required to reach equilibrium increased with decreasing 
a;*, ranging from ~ lOOte for a;* = 10~^ to ~ tB for a;* ^ 1, where tB = tb/cs is the 
Bondi time. In the equilibrium configuration, all the runs except cj* = built up dense 
material around the accretor. The material was not in a symmetric disk, but rather in a 
pinwheel shape. The velocity patterns associated with these pinwheels generally involved 
matter flowing in from the z direction through an accretion funnel and cither circulating or 
outflowing in the xy plane. This overall morphology agrees well with the results of Proga 
and Begelman (2003). Figures 3, 4, and 5 illustrate the morphology for lo^ — 10~^, 1, and 
lO^-^ 

A comparison of runs at different a;* values shows several clear trends. First, the accu- 
mulated pinwheel of matter is much smaller at large a;* than at small a;*. For a;* ^ 1, the 
accumulated matter extends out to ~ tb/S, while for a;* ^ 1 it is well inside the Bondi 
radius and is less massive. The size of the pinwheel changes only slightly with uj^ for a;* ^ 1, 
but shrinks dramatically for a;* ^ 1. Our theoretical description of the torus we expect for 
small a;* appears to be roughly correct: a torus of material accumulates that extends out to 
~ Tb and has a height comparable to its radius. The tori are chaotic and time-variable. For 
larger u;*, the pinwheels arc more regular and less chaotic, and are accompanied by leading 
and trailing shocks. The shocks are moderately strong, with Mach numbers ~ 2 — 3. 

For a;* ^ 1, we do not flnd that the torus has constant specific angular momentum, as 
Proga and Begelman (2003) did. In our torii, the specific angular momentum of the gas 
varies by orders of magnitude. It is unclear if this difference in results is due to different 
initial conditions, or due to the fact that our simulations are three-dimensional where theirs 
were two-dimensional. However, we agree with Proga and Begelman (2003) that the size and 
shape of the torus is largely independent of a;* (or the initial specific angular momentum /), 
and that its equilibrium structure arises due to the accumulation of material with too much 
vorticity (or I) to accrete. As Figures 3 and 4 show, accretion occurs only through a narrow 
funnel, and the (invariant) shape of the torus therefore determines the accretion rate. 
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4.2. 



Evolution of Circulation 



For each run we computed the average dimensionless circulation F as a function of time. 
We define the dimensionless circulation within a circle of radius re in the xy plane as 



where the line integral is evaluated on the circle r = r^. The circulation is a measure of 
the amount of vorticity within one Bondi radius of the accrctor. With this definition of 
dimensionless circulation, in the initial configuration of the velocity field the circulation is 
equal to the vorticity parameter, F* = a;*. Our prediction for the equilibrium circulation for 
a;* ^ 1 is F ~ 2'KCsrB, or F* ~ 2. 

We plot the circulation versus time along with the accretion rate versus time for each of 
our runs in Figure 6. The behavior of the circulation with time depends on the value of a;*. 
For a;* ^ 1, the circulation starts at a;* and gradually increases until it reaches ~ 0.4 — 0.8. 
Then it fluctuates around this mean. For a;* = 10°-^, F* decreases from its initial value 
until it again reaches the range ~ 0.4 — 0.8. As a;* increases, however, F* decreases less and 
less in the time it takes the accretion rate to reach equilibrium. The accretion rate and the 
circulation appear to be anti-correlated for lu^^ 1 in the period before equilibration: the 
initial buildup of circulation takes about the same amount of time as the initial decrease in 
the accretion rate. This supports our hypothesis that vortex lines are acting like magnetic 
flux lines and inhibiting accretion. 

For each run we estimate by eye the time at which the flow pattern and accretion rate 
reach equilibrium, and we compute the mean circulation after this time. We report the 
equilibration times and the mean circulations in Table 1 ; we plot the equilibrium value of F* 
in Figure 7, and the equilibration time in Figure 8. As the figures show, the equilibrium value 
of the circulation is at most very weakly dependent on the initial vorticity for a;* ^ 1, varying 
by less than a factor of 2 while the initial vorticity varies by a factor of 100. The flow pattern 
re-arranges itself to select F ~ 0.4 — 0.8 by accumulating a dense torus of material rotating 
at nearly Keplerian speeds around the accretor. The equilibrium value is smaller than our 
prediction of F* ~ 2, indicating that the rotational velocity is sub-Kcplcrian. This result is 
consistent with the morphology of the gas flow. The gas has a non-zero radial velocity and 
thus can be marginally bound even with tangential velocities that arc sub-Kcplcrian. Our 
analytic model, because it fails to predict the chaotic nature of the fiow pattern, overestimates 
the equilibrium value of F^ by a factor of ~ 4. It does, however, correctly predict that there 
is an equilibrium F* for a;* ^ 1. At a;* ^ 1, the pinwheel of unaccreted matter is confined 
closer to the accretor, and is unable to affect the circulation on length scales ~ tb. Thus, 
the circulation for larger cu^ stays constant. 




(35) 
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The equilibration time is also roughly consistent with our prediction of ^eq = 2iB/<^* 
for 0;=^ ^ 1. Part of the scatter of the line comes from the unavoidable subjectivity in our 
estimate of the equilibration time in such a chaotic flow. However, our model is clearly only 
approximate. While we correctly capture the trend of increasing equilibration time with 
decreasing u;*, our inability to model the turbulent flow means that our prediction of the 
equilibration time, like our prediction of F*, is accurate to at best a factor of a few. 

In the plots of and t^q versus (and as we show below, in the plot of M versus u^, as 
well), there seems to be a change in behavior between the cj* = 10^^ and u;* = 10^°'^ runs. 
Rather than smoothly interpolating bewteen c<j* ^ 1 and u^, ^ 1, the trend with uj^ appears 
to jump from one track to another. We cannot rule out the possibility that this is simply 
the result of chance and our sparse sampling a;* values. Near cu^ — 1, the equilibrium values 
of quantities like and t^q may be fluctuating chaotically with a;*, and the sharp jump we 
see from a;* = 10~^ to a;* = 10~^/^ may just be a fluctuation. However, it is also possible 
that there is a distinct regime of intermediate a;*, as well as the small and large a;*, and that 
several properties of the equilibrium flow change suddenly when one enters this new regime. 



4.3. Accretion Rates and Comparison to Theory 

We ran each of our simulations until the accretion rate converged to a steady state, and 
then we measured both the mean accretion rate and the fluctuations about the mean. Table 
1 summarizes our results and Figure 6 shows the accretion rate versus time for each run. We 
fllter out short time scale noise associated with the flnite resolution of our simulation by only 
computing the standard deviation in accretion rate averaged over 16 time steps. Nonetheless, 
discreteness error does lead to a non-zero standard deviation even in the case uj^, = 0, which 
is pure Bondi accretion and should have zero standard deviation. The standard deviation 
of 0.007% gives a lower hmit on our ability to measure fluctuations in the accretion rate. 
Since the standard deviations we measure for other runs are signiflcantly larger than this, 
they likely represent real physical time variability, not just simulation noise. 

We plot the mean accretion rate as a function of in Figure 9. We also show our 
theoretical prediction (solid line), which flts the data reasonably well. To obtain a somewhat 
better flt, we keep the function /(a;*) the same and do a least squares flt to obtain a pre- 
factor to replace 0.34 in equation (29). We flnd a best- flt value of 0.40. This value produces 
a good flt, as shown in the dashed line in Figure 9. The formula flts each of the simulation 
values to better than 40%. This is comparable to the accuracy of the Bondi-Hoyle formula 
at intermediate Mach numbers (Ruffert 1994; Ruffert and Arnett 1994). 
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4.4. Convergence Testing and Dependence on r^/r-Q 

To test the convergence of our simulations, we repeated one of our runs at higher 
resolution. Perna et al. (2003) have hypothesized that the accretion rate in Bondi accretion 
with vorticity is proportional to (rA/re)^, where p is between 0.5 and 1. Proga and Begelman 
(2003) saw a small dependence of the accretion rate on rx/r-Q for some of their simulations 
with smaller dynamic ranges, but the variation seemed to disappear in their simulations with 
larger dynamic ranges. In our sink particle method, the central object cannot accrete angular 
momentum, so in the absence of numerical viscosity, we would be in the limit ta/tb -C 1. 
However, as we have noted above, numerical viscosity does set an effective minimum size 
for our accretor. In the inner few cells of the simulation grid, where the circles in which the 
gas is attempting to move are poorly resolved, our code does not perfectly conserve angular 
momentum or vorticity. It is difficult to determine an equivalent accretor size set by this 
phenomenon, but we have found that even at our lowest resolution of re/ Ax = 40, our 
simulations produce disks and accumulation of circulation rather than Bondi-like flow for 
a;* as small as 10~^. Thus, the effective value of rx/r-Q set by our method must be ^ 10^^. 
(Recall that typical values of ta/tb in astrophysical systems are closer to 10~^.) Regardless 
of the true effective accretor radius imposed by numerical viscosity, by changing the grid size, 
and thus the rate of angular momentum dissipation, we are testing whether the accretion 
rate truly depends on ta/tb when ta/tb <C 1. 

To test for convergence, we re-ran our simulation with a;* = lO'''^ at a resolution of 
tb/Ax = 160, four times the fiducial resolution. All other aspects of the simulation were 
identical. We plot the accretion rate versus time for the two differently resolved versions of 
the run in Figure 10. As the results show, the exact shape of the accretion rate versus time 
graph is not the same at the different resolutions. This is to be expected in an unstable and 
chaotic flow. However, the overall accretion rates appear to be approximately the same. We 
find a mean accretion rate after equihbration of 0.088Mb for the run with tb/Ax = 40, and 
a mean of 0.11Mb for the run with tb/Ax = 160. The standard deviations in the accretion 
rates are 22% and 14%, respectively, while the difference in means is only 20%. Thus, our 
results are consistent with the hypothesis that there is no dependence of the accretion rate 
on ta/tb in the limit ta <^ tb- The results also suggest that our simulations are converged. 
In contrast, consider the hypothesis of Perna et al. (2003) that the accretion rate scales as 
{ta/tbY for P in the range 0.5 to 1. By increasing the resolution by a factor of 4, we should 
have decreased ta/tb by a factor of 4, and thus the Perna et al. (2003) proposal predicts 
we should have measured an accretion rate in the range 0.044Mb to 0.022Mb for the higher 
resolution simulation. This is clearly inconsistent with our results. We emphasize, however, 
that we have only tested the hydrodynamic case. Perna et al. (2003) made their hypothesis 
for both hydrodynamic and magnetohydrodynamic flows, and we have not tested the latter 
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case. 



5. Conclusions 

We provide a general theoretical framework for considering problems of accretion in 
a medium with vorticity. Using simple analytic estimates, we have been able to derive a 
formula for the rate of mass accretion onto a point particle of radius ta as a function of the 
vorticity present in the ambient medium. We define the dimensionless vorticity parameter 
a;* as 

uj^ = (36) 

Cs 

where tb is the accretor's Bondi radius and is the sound speed. We find that the accretion 
rate is ^ 

MM ~ 47rpoo^^^0.4/(a;,), (37) 

for ci;*^a;crit, where f{uj^) is the function defined by (27). We note that /(a;*) has the 
limiting behavior /(a;*) ~ 1 for a;* <^ 1, and 

/(a;,) ^ In (16a;,) (38) 
otto;* 

for uj^^ 1. Simulations show that our formulation provides a good fit to the overall shape of 
the curve of accretion rate versus vorticity. By calibrating our first-principle calculation using 
the results of our simulations, we provide a formula that agrees with the mean accretion rate 
we measure in our simulations to better than 40%, all over 3.5 orders of magnitude variation 
in vorticity. Our result is a natural extension of the Bondi-Hoyle-Lyttleton formula to 
rotating fiows. 

We are also able to roughly predict several other properties of the flow, including the 
overall morphology, the equilibrium vorticity, and the equilibration time. Our predictions for 
these quantities are considerably less accurate than for the accretion rate, but our formulae 
are correct to within factors of a few. We see ambiguous evidence for the existence of 
an intermediate uj^ regime, which is characterized by a rapid change in equilibrium flow 
properties starting between a;* = 0.1 and a;* = lO*^'^. Our analytic model does not predict 
such the existence of such a transition, and it is possible that what we have seen is a result 
of numerical noise. Regardless, it does seem clear that the analytic model runs into trouble 
when a;* ~ 1. Since we derived it for the regimes a;* <S 1 and a;* ^ 1, and simply interpolated 
between those, this is not surprising. 
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The simulations we use to check our formulation have very general and simple initial 
conditions. We have avoided complications arising from boundary conditions by moving the 
boundaries very far away from the accreting object, and our sink particle method allows 
us to avoid the need for an unrealistically large accretor radius. Our theoretical discussion 
shows that the accretion rate is controlled by a combination of vorticity build-up near the 
accretor and, for a;* ^ 1, the rate at which gas with low enough energy to accrete approaches 
the accreting particle. As a result, our results should be extremely robust. 

Our results allow direct application to observed systems. One may determine the vor- 
ticity in a astrophysical system by observing a velocity gradient across it, which is often 
feasible using radio line observations. Consider an example in the context of star formation: 
Goodman et al (1993) observed dense molecular cores in nearby star-forming regions. They 
found velocity gradients of ~ 0.3 — 3 km s~^ pc~^, in regions with typical temperatures of 
~ 10 K (sound speed Cg ~ 0.2 km s~^ for a mean molecular mass of 2.33mH). This means 
that a one solar mass protostar (re ~ 0.1 pc) inside such a dense core has a vorticity pa- 
rameter LV^, ~ 0.15 — 1.5 (specific angular momentum / ~ 8 x 10^0 — 8 x 10^^ gr cm s~^, and 
should have its accretion rate reduced by a factor of ~ 3 — 7 relative to the Bondi rate. This 
result has potential implications for models of star formation in which protostars gain mass 
through a process of competitive accretion. 
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A. Approximation of /(a;*) for Large a;* 



We wish to derive an approximate formula for 




(Al) 
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in the limit a;* — > oo. We begin by breaking it into two terms, 



1 



At the upper hmit of the integration, y' — y'^g^i the two terms under the square root are 
equal and the integrand vanishes. For smaller values of y' , the first term goes to 4, while the 
second term scales as y'^. Thus, the second term is significant only near the upper limit of 
integration. Examining Figure 1, we can see that the region y' ^ l/max niakes only a small 
contribution to the integral in the large uj^ case. Thus, we may obtain an approximation for 
/(a;*) when a;* 3> 1 by dropping the second term, giving 



1 /■!/max{'^*) / 4 

2 

■smh-'[uj,y'^^{uj.)]. (A4) 



The upper limit of integration is 

2/3 



vLa^M = : . 1/3 — , (A5) 



3V3 (gu;, + ,/3T8h4) - 32/3 
which for a;* ^ 1 behaves as 

/ 2 \ 

lim y'^^M = — . (A6) 

Substituting into (A4), we find 

fM ^ sinh-i f(2u;*)'^'l ■ (A7) 
TTo;* L J 

Taking the limit a;* — > cxd, this becomes 

fM ~ In (16a;*) . (A8) 
otto;* 

Numerical integration shows that this approximation is fairly accurate over a wide range of 
a;*, but the convergence is logarithmic and the approximation does not become very good 
until uj^ is quite large. It is within 10% for cj* > 4 and within 5% for a;* > 53, but does not 
become accurate to 1% until a;* > 1.8 x 10^^. 
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Fig. 1. — The plot shows the area A{uj^), as defined by equation (21), for different values of 
a;*. The numbers indicate the value of logo;* for each curve. 
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Fig. 2. — The solid shows f{uj*), as defined by equation (27). The dashed hne shows the 
large a;* approximation given by equation (28). 
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Table 1. Simulation parameters and results. 



a;* 


AXmin/rB 


grange /'"B 


^run / 


teq/tB 


M/Mb 









0.025 


[-40, 40] 


77 


7 


0.99 


0.007 


2.3 X 10-^ 


10-2 


0.025 


[-40, 40] 


200 


100 


0.36 


0.099 


0.49 


10-1-5 


0.025 


[-40, 40] 


120 


30 


0.31 


0.099 


0.56 


10-1 


0.025 


[-40, 40] 


50 


15 


0.25 


0.061 


0.65 


10-0.5 


0.025 


[-40, 40] 


63 


30 


0.34 


0.15 


0.38 


10° 


0.025 


[-40, 40] 


19 


12 


0.22 


0.10 


0.42 


100.5 


0.025 


[-12.5,12.5] 


18 


4 


0.088 


0.22 


0.75 




0.0061 


[-6.3,6.3] 


4.3 


1.5 


0.048 


0.19 


6.7 


101.5 


0.0061 


[-12.5,12.5] 


0.96 


0.7 


0.018 


0.037 


29.2 



Note. — Col. (1): Vorticity parameter. Col. (2): Grid spacing in Bondi radii. 
Col. (3): Size of simulation region in the y direction. Col. (4): Run duration in 
Bondi times. Col. (5): Approximate time the accretion rate reaches equilibrium, in 
Bondi times. Col. (6): Mean accretion rate, in units of the Bondi rate. Col. (7): 
Standard deviation in accretion rate, as a fraction of the mean accretion rate. Col. 
(8): Mean dimensionless circulation. 
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Fig. 3. — The plots show shces through the origin in the xy and xz planes for the tu* = 
10^^ run. The grayscale indicates log density, and the arrows indicate velocity. Note the 
differences in scale from Figures 4 and 5. 



- 28 - 




2.50 



2.08 



1.67 



1.25 



0.83 



0.42 



I 0.00 
1.0 10 00 Log (p/pj 




1.0 



0.5 



0.0 



-0.5 



1.0 



-~-^^^^-v^-^NSS\\\\\\\\ \ \ I I I 
^^^^--•^NNSS\\\\\\\,'^ \ \ V i I 

~^^^v.^-,^...^s.sSN\\\S,S,\,\, \ \ I I \ 

- ^ ^ ^^^^^^-^•^-^•^•^•^WWW \ ^ \ 1 I 

\ ^^^NS\\NNN.W\\\\\\\\\ \ \ \ 
' ' ' ' ' ' ^ "'^^\\\\\\\\\\\\\\\\ 

i I V \ wwwwwwww 

.www I t t 1 1 f ? f^JJJJ t t t 

.WW 1 I / /////■///////////// 
t I I /////////////////////// 

^^^^^^^^^y^/'y/// / / / / / / / t t t 

^^^^^^^^//////fff 1 1 1 1 t I 

^^^^^^^^^^/■/■/////^rtit t I t I 
^^^^^^^^^^//■/'////f/ttti t t t 
^^^^^^^^^.'.'^////fftit t r t I I 



\\\\\\\\\\\\ 
\ \ \ \ \ \ www 

\W SWWNSSV 
\ \ S \N\\NS.NS.N 



/ / X /■/■,'/' y y , 
! t t r / / y / / / / 

t t T t t r ; / / / / 
t f t t t w t f f ; 



1.0 



-0.5 



0.0 

X / R 



0.5 



B 



1.0 10.00 



2.50 



2.08 



1 .67 



1.25 



0.83 



0.42 




0.00 

Log {p/pj 



Fig. 4. — The plots show shces through the origin in the xy and xz planes for the u^, = 1 run. 
The grayscale indicates log density, and the arrows indicate velocity. Note the differences in 
scale from Figures 3 and 5. 
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Fig. 5. — The plots show shces through the origin in the xy and xz planes for the tu* = 
10^'^ run. The grayscale indicates log density, and the arrows indicate velocity. Note the 
differences in scale from Figures 3 and 4. 
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Fig. 6. — The plots show the accretion rate versus time (sohd hnes) and circulation versus 
time (dashed lines) for the run with the value of uj^ indicated in the panel. The dotted 
vertical lines indicate the time at which we conclude the system has reached equilibrium. 
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Fig. 7. — The asterisks indicate the mean value of after equihbration for each of our runs 
as a function of a;*. 
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Fig. 8. — The asterisks indicate the time (in units of t^) at which the run reaches equihbrium 
for each value of a;*. The hne is our prediction tgq = Ste/*^*- 
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Fig. 9. — The plot shows mean accretion rate, in units of the Bondi rate, versus uj^. The 
asterisks are the simulation data, with error bars indicating two times the standard deviation 
in the accretion rate. The solid line is our theoretical prediction (29). The dashed hne is 
the theoretical prediction scaled by a constant factor to give the best possible fit to the 
simulation data. 
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Fig. 10. — The plot shows the accretion rate versus time for our convergence test. The thick 
hne is the a;* = 10°"^ run at a resolution of tb/Ax = 40, and the thin line is the run at a 
resolution of tb/Ax = 160. 



